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Cooled, low-loss nanomechanical resonators offer the prospect of directly observing the quantum 
dynamics of mesoscopic systems. However, the present state of the art requires cooling down to the 
milliKelvin regime in order to observe quantum effects. Here we present an active feedback strategy 
based on continuous observation of the resonator position for the purpose of obtaining these low 
temperatures. In addition, we apply this to an experimentally realizable configuration, where the 
position monitoring is carried out by a single-electron transistor. Our estimates indicate that with 
current technology this technique is likely to bring the required low temperatures within reach. 



PACS numbers: 85.85.+j,85.35.Gv,03.65.Ta,45.80.+r 



I. INTRODUCTION 



Nanomechanical resonators are now being built with 
quality factors in the range, Q f=a 10 , and resonance fre- 
quencies of up to several hundred MHai. The ground 
state energy of these devices can correspond to tempera- 
tures in the milliKelvin range. As a result, the observa- 
tion of quantum behavior in these devices is becoming a 
real possibility 2 -. To detect such behavior, the resonator 
must be sufficiently cold; since a quantum harmonic os- 
cillator driven by thermal noise behaves as a classical 
oscillator driven by thermal noise, one must ensure that 
the signatures of quantum effects are not swamped by the 
thermal behavior. The approach taken so far to achieve 
low temperatures is to place the resonator in a refriger- 
ator. However, cooling very small devices in this way is 
inherently inefficient in that the system becomes weakly 
coupled to the thermal bath. Here we explore the possi- 
bility of using feedback control to effect 'active' cooling 
of the resonator, in order to cool below the possible limits 
set by the 'passive' refrigeration technique. 

To perform such feedback cooling the resonator must 
be monitored, and the result fed back in real time to af- 
fect the dynamics. A practical method of performing a 
continuous measurement of the position of the resonator 
is to use a Single-Electron Transistor (SET)iM£. To mea- 
sure the position of the resonator one locates the central 
island of the SET next to the resonator. When the res- 
onator is charged, and the SET is biased so that current 
flows through it, changes in the resonator's position alter 
the potential on the central island, which in turn changes 
the current. The current therefore provides a continuous 
measurement of the position of the resonator, and this is 
just what is required for implementing a linear feedback 
cooling algorithm^. A feedback force can be applied by 
applying a voltage to a gate capacitively coupled to the 
resonator, and adjusting the voltage so as to damp the 
resonator (see Fig. 1), or by passing a variable current 
through the oscillator in the presence of a fixed external 
magnetic field. We will analyze the first system, although 
the results should apply to the second as well. In our 



analysis we will use the theory of the dc-SET. While an 
experiment would most likely use a radio-frequency SET, 
the characteristic frequency of a SET is typically of the 
order of 10 GHz, so that the RF drive looks constant to 
the SET, and the dc-SET equations can be used. 

We will use a quantum mechanical model of the mea- 
surement and feedback process, but discuss how in this 
case, such a description is equivalent to a classical mea- 
surement of a noisy classical system. Thus, this article 
is intended for both experimentalists familiar with clas- 
sical descriptions of noise in systems, as well as quantum 
measurement theorists. 

Rather than performing a microscopic analysis of the 
measurement process in terms of the interaction of the 
SET and the resonator, we start by introducing equations 
which describe the continuous observation of a quantum 
observable, and show how this includes the shot noise 
and back-action, these being the key sources of noise 
in a continuous quantum measurement. This descrip- 
tion can then be tailored to the case of a measurement 
with a SET by choosing the parameters so that the noise 
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FIG. 1: A schematic of the resonator, measuring and feed- 
back apparatus. As the resonator moves closer to the SET, 
the current flowing through the SET changes, and that infor- 
mation is then used to generate a feedback voltage applied to 
an actuating gate. 
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sources match those calculated in microscopic noise anal- 
yses which have been performed for the SET 3 - 8 . 

A treatment of the continuous quantum measurement 
of a two-state system using a SET has been carried out 
by Korotkov^, using what might be referred to as a par- 
tially microscopic approach. The equations we use here 
may be derived by replacing the two-state observable in 
those equations by the resonator position. A full analysis, 
along the lines of those performed for quantum optical 
system s 10 ! 11 , can also be expected to produce the same 
equations under reasonable approximations. The form 
of these equations is determined by how information is 
obtained, and not by the specific implementation, which 
explains why the form of the equations is similar in op- 
tical position measurements and position measurement 
using SETs. If the measurement is of a physical observ- 
able, and the resulting error about the expectation value 
of that observable in a short time interval At is Gaus- 
sian, then the most straightforward implementation of 
that measurement process has the form used here. 

In the next section we introduce the equations that 
describe a continuous measurement process, derive the 
form of the resulting noise, and give the equivalent clas- 
sical model. We then discuss how this model can be 
applied to position measurement using a SET, and com- 
pare our formulae to those derived using a semi-classical 
treatment of the SET 3 ^ in order to express our results 
in terms of experimental parameters. In Section III we 
discuss the implementation of a feedback algorithm and 
calculate the minimum achievable temperature in terms 
of physical parameters. We then calculate estimates of 
realistic achievable temperatures for an an experimen- 
tally realizable sample system in Section IV, and finally 
conclude with a summary of the results obtained. 



II. CONTINUOUS QUANTUM 
MEASUREMENT OF POSITION 

Given a quantum system whose state is specified by 
the density matrix p, and whose evolution is determined 
by the Hamiltonian H, then a continuous measurement 
of the observable O of that system, which provides the 
continuous output results (measurement record), 



dr = (0)dt 
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dW, 



(1) 



induces the following evolution of the system&iSii 3 . 

dp = -i[H,p]dt-k[0,[0,p]]dt 

+ V2k(Op + pO -2(0))dW. (2) 

Here k is proportional to the measurement strength and 
dW is a Weiner process. The noise contained in the mea- 
surement record is a necessary result of the fact that only 
a finite amount of information is obtained regarding the 
observable O in a finite time. This direct noise on the 
record is called the shot noise. However, this is not the 



only noise resulting from the measurement process. As 
a result of Heisenberg's uncertainty relation, information 
about one observable makes other observables less cer- 
tain. Due to the dynamics, the uncertainty (noise) in 
these observables can feed into the observable being mea- 
sured. This source of noise is referred to as back-action. 
If the Hamiltonian is such that the increased uncertainty 
is not fed back into the observable being measured, then 
the measurement is referred to as 'back-action evading.' 

Now let us examine the case of a position measurement 
on a harmonic oscillator. To do this, we set O = x, and 
the Hamiltonian becomes 



H 



P 

2m 



1 2 2 

-muj x , 



(3) 



where m is the mass of the particle, luq is the (angular) 
frequency of the oscillation, and x and p are the posi- 
tion and momentum operators, respectively. To make our 
model sufficiently realistic, we need to include two more 
sources of noise: the first is the intrinsic thermal noise of 
the harmonic oscillator, and the second is the possibility 
that the oscillator may be driven by white noise over and 
above that required by Heisenberg's uncertainty principle 
(excess 'technical noise'). 

The second of these is easily included by adding a 
term — 0[x, [x,p\] to the equation of motion of p; this 
describes a noise term identical to the one caused by the 
back-action, but without the corresponding dynamics of 
p associated with obtaining a measurement result which 
causes the back-action. It is equivalent to adding a term 
linear in x to the Hamiltonian J2J multiplied by white 
noise. 

The inclusion of thermal fluctuations is only a little 
more involved, and can be achieved by coupling the os- 
cillator to a thermal bath. In our case the effect of the 
thermal bath may be included by adding the 'standard 
Brownian motion master equation' (SBMME) 14 to our 
equation of motion for p: 



dp 



--[H,p]dt- — [x,{p,p}+]dt 
- (k + ^+^fcoth^j [ X ,[x,p]]dt 



27; 



2k(xp + px-2(x)p)dW 



(4) 



where r = ^, Q being the quality factor of the res- 
onator. The two terms proportional to T are due to the 
inclusion of the SBMME, the first representing dissipa- 
tion due to the reservoir while the second is a diffusion 
term due to environmental fluctuations. Here we are us- 
ing an approximate form of the SBMME appropriate for 
the weak coupling regime (small T, large Q) but covering 
all ranges of temperatures^. Since the nanomechanical 
resonators we consider all have large values of Q, the 
weak coupling requirement is easily satisfied. The tem- 
perature dependence of the diffusion coefficient is given 
by coth(7ia;o/2fcBT) so that the diffusion does not van- 
ish as ksT — > 0: this correctly accounts for the exis- 
tence of quantum vacuum fluctuations which exist even 
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at zero temperature. In the absence of a rigorous charac- 
terization of the dissipation channels of nanomechanical 
systems there is as yet no need to include a more sophis- 
ticated description of SBMME environmental effects^. 
Phcnomenological corrections to the SBMME such as 
temperature dependence of T can be added if needed, 
but these are not significant effects in the high-Q regime. 

We also need to include in our model the possibility 
that there is noise driving the oscillator which is cor- 
related with the noise on the measurement record (the 
shot noise). This can happen if the noisy behavior of 
the oscillator explicitly causes some of the noise in the 
measurement apparatus, or vice versa. In this situation, 
the measurement record contains more information about 
the oscillator position, so when it comes to adding feed- 
back, we are able to cool the oscillator further than would 
otherwise be expected. In Eq. J5J the noise driving the 
oscillator is purely the quantum back-action. It may ap- 
pear from Eqs. and J3J) that the quantum back-action 
is correlated with the shot noise due to the fact that 
the same noise term (dW) appears in both equations. 
However, this is not the case. The term proportional to 
dW which appears in the equation for p describes the 
random way in which the measurement changes the ob- 
servers state of knowledge about the system. Thus, in 
general, this noise term decreases the entropy of p. The 
back-action noise, which is driving the oscillator and con- 
sequently increasing the entropy of p, is described by the 
term proportional to k. The quantum back-action is, in 
fact, completely uncorrelated with the shot noise. 

To drive the oscillator with a random force, one applies 
the Hamiltonian ht;(t)x, where £(t) is the magnitude of 
the random force. We can choose to be correlated 
with the shot noise, with correlation coefficient «, by set- 
ting 



d£ = V2^(^dW + VT^ndV), 



(5) 



where dVis a Wiener noise uncorrelated with dW . The 
resulting spectral density of £(t) is a, so that (£(£)£(£')) = 
aS(t—t'). The Stratonovich equation which describes the 
driving by is 

|$ = -*e(t)a#), (6) 
and converting this to an Ito equation gives 

d\i>) = -iV2ax\ip)d£ - ax 2 \ip)dt. (7) 

Converting the Ito equation further to an equation for p 
one obtains 

(8) 



dp = —a[x, [x, p]]dt — i\J2a[x, p]d^ 



Since the observer has access to dW, but not to dV, she 
must average over dV, and this gives 



dp = —a[x, [x, p]]dt — iy2na[x, p]dW. 



(9) 



If we allow part of the excess noise given by [3 in our 
model to be due to driving by the shot noise dW (that 



is, this noise is correlated with the shot noise dW with 
correlation coefficient k) then the equation of motion for 
the system becomes 



dp 



■-[H,p]dt- — [x,{p,p} + ]dt 



k + fj + 



■mojoT hu>o 
coth 



2k B T 



2% 

- iy/2K0[x, p]dW 

+ V2k(xp + px - 2(x)p)dW 



[x, [x,p]]dt 



(10) 



This completes our quantum mechanical description of a 
resonator under continuous observation. 

Now that we have an equation that includes all the 
relevant noise terms, the noise spectrum of the measure- 
ment record can be obtained: 
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8k 



k + (3 + 
2{h/m) 2 



mu T hu \ 

2h °° 2k B T J 



r 2 < 



J 2 + (ij 2 - U! 2 ) 2 



(11) 



The first term is the shot noise, which is white, the term 
proportional to k is the quantum back-action, the term 
proportional to T is the effect of the noise from the res- 
onator's thermal environment, and the term proportional 
to (3 gives any excess noise over and above the necessary 
quantum back-action. Note that the last three terms all 
have the same form as a function of oj. This is because 
they are all white noises filtered through the harmonic 
oscillator spectral response function. 

While our treatment so far has been fully quantum 
mechanical, it is worth noting that a completely classical 
model of a measured, damped oscillator will reproduce 
the dynamics of this measured quantum system, so long 
as the initial density matrix is Gaussian in x and j&. 
Thus, one can understand the behavior of the oscillator 
in terms of classical noise and a classical measurement 
process. The equations of motion for the position, x c , 
and momentum p c of this equivalent classical oscillator 



dx c 
dp c 



= —p c dt 



(12) 



= -mu$x c dt - Tp c dt + hV2kdY c + hy/2f3dV c 



mhujnT coth —dU r 

2k B T 



(13) 



where dY c , dV c and dll c are each zero- mean Gaussian 
white noise, and mutually uncorrelated. The position of 
the oscillator is then observed by a continuous classical 
measurement, which generates the output record 



dr c = x c dt + —dZ, 



(14) 



and where dZ c is zero-mean Gaussian white noise, uncor- 
related with dY c . The noise term dY c is what is required 
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in the classical model to correctly include the back-action 
of the quantum measurement process. It is now explicit 
that this noise is uncorrelated with the shot noise on the 
measurement, dZ c . 

In the classical case, the observer's state of knowledge 
about the oscillator is given by a joint probability density 
over x c and p c - This probability density is the classical 
equivalent of the density matrix p. So long as the initial 
probability density is Gaussian, it remains Gaussian as 
time passes, and as a result the observer's full state of 
knowledge may be represented by merely 5 variables: the 
mean position and momentum, (x c ) and (p c ), and the 
variances and covariance, given by 



(xl) - (x c )\ 



Xp 



{XcPc 



- (Xc)(Pc 



(15) 
(16) 
(17) 



It is the means (x c ) and (p c ) (being the observer's best 
estimates of the value of x c and p c ) which are the clas- 
sical equivalents of the quantum expectation values (a;) 
and (p) . It turns out that if one writes the classical mea- 
surement record as 



dr c = (x c )dt + ^-rdW c , 
8k 



(18) 



then dW c is zero-mean Gaussian white noisei, uncorre- 
lated with dZ c . The classical model is then equivalent to 
the quantum model if we equate dW c with the quantum 
measurement noise, dW, and correlate dV c with dW c , so 
that (V c (t)dW c {t')) = n8{t-t'). 



III. CONTINUOUS MEASUREMENT WITH A 
SINGLE ELECTRON TRANSISTOR 

Having obtained a model which is sufficiently general 
to encompass the dynamics of a resonator monitored by 
a SET, we need to express the theoretical parameters k, 
(3, and k in terms of the actual experimental parameters 
of the SET. Since it is by measuring current through 
the SET that we measure the resonator position, it is 
the spectral density of this current which determines the 
shot noise of the measurement. The back-action from 
the measurement is due to the action of the SET on the 
resonator, which is the force that the resonator feels from 
the charge on the SET island. As a result the back-action 
noise, /3, can be calculated from the spectral density of 
the charge fluctuations on the SET island, and hence k 
is determined by the correlation between the current and 
the island charge fluctuations. 

However, the dynamics of the SET are sufficiently 
complex that analytic results for these spectra have as 
yet only been obtained for certain parameter regimes. 
These calculations have been performed by Zhang and 
BlencoweSi, using previous results of Korotko\^. The 
technique used is to approximate the dynamics of the 



electron tunneling on and off the SET island by a classi- 
cal master equation. That is, the electrons are assumed 
to tunnel independently across each of the junctions, with 
certain rates (the rates being obtained using a perturba- 
tive quantum calculation). This ignores the possibility 
that electrons will tunnel coherently across both junc- 
tions simultaneously, a quantum effect referred to as co- 
tunneling. 

However, it is important to note that the above "classi- 
cal" method for calculating the charge fluctuations does 
not include the quantum back-action noise. This can 
be seen from the following argument— ■. In the classical 
treatment, since the fluctuating force on the resonator is 
due to the electrons jumping on and off the island, in 
principle the time history of this force can be known by 
detecting the electrons flowing in the circuit. In prin- 
ciple then, the effect of the noise can be known, and if 
desired, undone. As a result it cannot include the quan- 
tum back-action, since this cannot, even in principle, be 
undone. 

The approach we will take here is to use the classical 
calculation of the force noise on the resonator, which de- 
termines /?, and add to it the necessary back-action as 
required by quantum mechanics, which is given by k. In 
doing this we note that the ratio k/(3 provides a diag- 
nostic tool for determining when the classical calculation 
breaks down; if k/f3 -C 1 is not satisfied, then the clas- 
sical calculation no longer provides a good estimate of 
the total force noise on the resonator. Thus it should 
be noted that if k/(3 > 1, then the classical calculation 
cannot be relied upon. That is, it is possible in this case 
that the total noise on the resonator is significantly larger 
than our estimate k + [3, due to quantum contributions 
not taken into account in the classical calculation. 

We find that in the regions of best cooling, which we 
explore in the following, fc is larger than (3 (although 
near-optimal cooling can be obtained with k < (3, and in 
particular we will give as an example results for k = (3/ A). 
Hence our calculations should be regarded as estimates of 
the performance of the feedback algorithm, rather than 
exact results. We note, however, that a more sophisti- 
cated analysis using the diagrammatic techniques devel- 
oped by Schoeller and Schonii might provide analytic, or 
semi-analytic results for the parameter regime of most in- 
terest for quantum measurement and control, and there- 
fore may provide a method for more accurate calcula- 
tions. 

The spectral densities given by the classical calculation 
are derived in Appendix ^ Approximations which are 
used in the derivation are detailed there, and come pri- 
marily from Zhang and Blencowe^. The noise spectrum 
of the displacement of the resonator due to the shot noise 
of the SET current is 



Si 



5i(w) 



x 



(dlds/dx) 



(19) 



where Si(uj) is the spectral density of the shot noise, 
given in Eq. (|A10|) . and / is the current through the SET, 
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given in Eq. I|A8(1 . The dependence of the current on the 
displacement of the resonator comes from its dependence 
on the gate capacitance, which can be approximated by 



(20) 



The shot noise, Si{uj) is, to a very good approximation, 
frequency independent, as required by our quantum mea- 
surement model. Thus 



8k~ xlu=0 ~ (dl ds /dx) 2 



(21) 



The spectral density of the classical part of the displace- 
ment noise due to the fluctuating force on the resonator 
is 



Sx(o>) 



S F (to)/m 2 



T 2 u 2 + (lu 2 - uj 2 ) 



2\2 ' 



(22) 



where Sf{w) is the spectral density of the fluctuating 
force given in Eq. l|A13Jl . Since, once again, Sf{w) is 
effectively frequency independent, we have 



2h 2 



(23) 



The correlation coefficient, k, between the shot noise and 
the excess back action is therefore simply the correlation, 
C, between Si and Sp, which is given in Eq. JMTJl. 



IV. FEEDBACK CONTROL 

We wish to cool the dynamics of the resonator by using 
the information obtained continuously about the state of 
the resonator to direct a time-dependent external force. 
Such a force may be applied, for example, by passing cur- 
rent through the resonator and immersing it in a mag- 
netic field. It can also be applied by placing an actuating 
gate near the resonator, and varying the potential dif- 
ference between the charged resonator and the actuating 
gate. 

In this case the results of modern optimal control the- 
ory apply, since the dynamics of the resonator are equiv- 
alent to that of a classical oscillator driven by Gaussian 
noise, so long as we restrict ourselves to a linear external 
forc o 6 ' 20 . This allows us to obtain the optimal feedback 
algorithm in a straightforward manner. Choosing the 
minimization of the energy of the resonator as the feed- 
back objective it turns out that as long as the force we 
apply is sufficiently large, this force should be chosen to 



To calculate the average energy of the controlled res- 
onator, we first need the equations of motion for the mean 
values in the continually observed and controlled case, 
which arc 



d(x) = —dt~T{x)dt + 2V2ka 2 x dW 



(25) 



d(p) = —muj 2 (x)dt — 2T(p)dt — j(mui(x) + (p))dt 



'2KpndW + 2\/2kai v dW 1 



and for the covariances, 
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(26) 

(27) 
(28) 

(29) 



In these equations, a 2 and a 2 are the variances in position 
and momentum, respectively, and 



G x V = ^(xp + px) - (x)(p) 



(30) 



is the symmetrized covariance. This system of equations 
is exactly equivalent to Eq. (|10|l as long as the initial 
state is Gaussian. In order to solve this set of equations 
most easily, we make what we call the truncated Gaus- 
sian approximation. We assume that the feedback rate 
7 is much larger than the system's small intrinsic damp- 
ing r, and we therefore drop all terms proportional to T 
from the above equations. This approximation is easily 
justified for current experiments. 

The steady state solutions to these equations are 



V2u 
~8k 



V2m^u} 



A. 



2, ,3 



+ 



8k 

muj 



A + A 3/2 ] 
y/ nfikVA, 



8k 



-A. 



(31) 

(32) 
(33) 



where 



F = -j(muJo(x) + (p)), 



(24) 



where 7 is a rate constant which determines the overall 
strength of the force. This equation gives optimal per- 
formance so long as 7 3> which is within reach of 
current experiments, as detailed below. 



A + l = 



1 + 16- 



kh 2 [k + (1 - n)p + m&r coth ^L)} 



1/2 
(34) 



In the limit of both large and small values of k, A ~ k. 
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Using the steady state solutions for the variances a 2 
and a 2 p in Eqs. it25)l and lf^|) . and calculating the vari- 
ances of (x) and (p) , we obtain 
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(35) 



(36) 



The average steady-state energy of the oscillator un- 
der observation and feedback is the sum of the intrinsic 
variances of the Gaussian steady state and the variances 
in the observer's measurement record. Thus 
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1 



A I A 1 A \ 1 

-muj (a x + a {x) ) + 



2m 



V2A + A + — A 3 / 2 + — A 2 
2 47 
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Here we have used the simplifying assumption 7 ^> u>, 
since this is inherent in the optimal control condition. 

It is clear from Eqs. Ij34(l and (|38J) that reducing the 
background temperature allows for lower final tempera- 
tures. Extremely low values of k lead to heating as can be 
seen from the fact that A ~ k. For large k (correspond- 
ing to large gate voltage), the increased sensitivity of the 
measurement cancels the increased disturbance due to 
the measurement, with the result that the minimal tem- 
perature levels off as k is increased. 



V. ESTIMATES FOR ACHIEVABLE 
TEMPERATURES 

Current refrigeration technology allows experiments on 
nanomechanical resonators to be performed at tempera- 
tures of about 100 mK. It is therefore sensible to as- 
sume that the feedback algorithm will be applied to a 
device which is initially at this temperature. In such ex- 
periments the resonators typically have fundamental fre- 
quencies in the range /o = 1 — 100 MHz. As our example 
system we take a realistic resonator with fo = 12 MHz, 
which is 6 /zm in length, 50 nm wide, and 150 nm thick. 
We restrict ourselves to relatively low frequencies because 
of the limits of feedback circuitry, which we estimate can 
easily operate at 50 MHz. The effective mass of such 
a resonator is roughly 10 -16 kg. An achievable quality 
factor, Q, is on the order of 10 4 . 
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FIG. 2: The steady-state average occupation number, (A), as 
a function of the gate voltage (solid line), plotted along with 
the ratio k//3 (dashed line), and the drain-source current, 
(dot-dash line). The lower dotted line gives the minima of 
(N), and the upper dotted line gives the value of {N) when 
fe//3 = l/4. 



Realistic values for the resistances and capacitances of 
the junctions of a SET which would be used to monitor 
the resonator are R\ = Ri = 50 kil, and C\ — C2 = 
100 aF, and we place it d ~ 100 nm from the resonator. 
We estimate that the capacitance between the gate of the 
SET and the resonator will be roughly C g = 50 aF, so 
that C s = 250 aF (C s = 2C 3 + C g ). It is° important to 
note that the analysis we use in the appendix to obtain 
the noise spectra is only a good approximation in certain 
parameter regimes. In particular, we require that Vd s , 
being the drain-source voltage across the SET, satisfies 
U ds < e/C s . 

To apply the feedback force, we place the resonator 
100 nm from the actuating gate, and allow the controller 
to vary the voltage difference between the gate and the 
resonator between —4 V to 4 V. The capacitance of this 
arrangement is about 50 aF, so the maximum force that 
can be applied to the resonator is of the order of 10~ 8 N. 
This corresponds to 7 rj 1.08 x 10 13 s _1 , which is much 
larger than to and T, as required by the optimal control 
condition and truncated Gaussian approximation used in 
the previous section. 

Before giving theoretical estimates for the achievable 
steady-state temperature (or equivalently, the steady- 
state average occupation number of the oscillator, (N) = 
(a* a)), we need to explain two subtleties which affect the 
presentation of our results. When one examines the de- 
pendence of the steady-state (N) on the gate voltage, 
once finds that it oscillates very rapidly, with minima 
occurring in closely spaced pairs. Since V g is experimen- 
tally easy to tune, all else being equal it would make 
sense simply to plot these minima and ignore the com- 
plex structure. However, the locations of the minima are 
such that k/(3 can be large, and as discussed in section lTTll 
our results are more trustworthy the smaller k/f3. The 
situation is shown in detail in Figure [3 in which we dis- 
play, as a function of Vg, two pairs of the (N) minima, 
as well as k/ f3 and the current Ids- When presenting re- 
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FIG. 3: Estimates for minimum achievable temperatures as 
a function of gate voltage for a range of initial temperatures. 
The dotted lines give the minimum temperature under the 
additional restriction that /3 = 4k. From top to bottom, the 
initial temperatures are 2 K, 1 K, 500 mK, and 100 mK. 



suits in what follows, we will plot both the minima with 
respect to V g , and the temperature which results if we 
demand that k//3 < 1 /4. For clarity these points are also 
displayed in Figure [3 

As an example of the relative magnitudes of the various 
noise sources at the minima displayed in Figure [21 if we 
set the drain-source voltage at Vd s = e/(4Cs) = 0.16 mV 
and the gate voltage at V g ~ 1 V, then the noise sources 
are: 



-,29 ™-2 -l 



/3 = 4.47 x ltr- m S 



k = 4.2/3 
mu>oT 



2h 



coth —4= = 209 13 
2k B T 



(39) 
(40) 

(41) 



and the correlation coefficient is k = 0.633. 

Using the above parameter values to calculate the 
minimum achievable temperature, we find that A = 
4.95 x 10~ 5 , and T = 2.06 mK. This corresponds to 
an energy of about E ss = 2.85 x 10~ 26 J, and an aver- 
age occupation number (N) = 3.09. While this is very 
encouraging, ideally one wants to cool below the energy 
of the first excited state, and we now examine what is 
required to do this. 

While classically an increase in measurement strength 
would automatically lead to improved tracking of the 
resonator, and therefore more efficient cooling, quantum 
mechanically the situation is more complex due to the 
fact that more precise measurement also leads to in- 
creased heating due to back-action. Nevertheless, In the 
present case one finds that the increased sensitivity of 
the measurement with increasing measurement strength 
effectively cancels this heating, and as a result larger V g 
corresponds to better cooling. However, after a sharp 




FIG. 4: Estimates for minimum achievable temperatures as a 
function of gate voltage for a range of resonator quality factors 
and an initial temperature of 100 mK. The dotted lines give 
the minimum temperature under the additional restriction 
that P = 4fc. From top to bottom, the quality factors are 10 3 , 
10 4 , 10 s , and 10 6 . A quality factor of 10 4 is achievable with 
current technology. 



increase in cooling with increasing V s , the minimal tem- 
perature levels off, so greater V g no longer provides much 
benefit. In addition, at some V g snap-in is likely to oc- 
cur as the force between the SET gate and the resonator 
becomes too strong. This voltage, in our example sys- 
tem, is estimated to be roughly 4 V. As a result, we 
limit ourselves to V g < 4 V. At V g = 4 V the steady 
state energy E — 8.37 x 10 - 27 J, which is below the 
energy of the first excited state. This corresponds to 
T = 0.61 mK, and (N) = 0.555. Thus, if the energy 
were to be measured directly, immediately after turning 
off the feedback, energy jumps as a signature of quan- 
tum behavior may well be observable. As an indication 
of the return from increasing the gate voltage, the steady 
state energy is E = 1.48 X 10 -26 J for V g ps 2 V, which 
corresponds to (N) — 1.36. 

In Figure [3] we plot the theoretical estimates for the 
achievable steady-state temperature as a function of V g 
for a range of starting temperatures. The solid lines cor- 
respond to the absolute minima, and the dotted lines to 
the minima under the restriction that k < /3/4. Of par- 
ticular interest is the fact that for a starting temperature 
of 2 K (ie. with pumped liquid He), we obtain minimum 
temperatures in the range of 50 mK. Thus, even for an 
initial temperature of 2 K, feedback cooling might well 
be able to compete with dilution refrigerators. If the res- 
onator is first cooled in a dilution refrigerator, and then 
feedback cooled, then the semi-classical theory predicts 
achievable temperatures below 1 mK. In Figure 0] we 
plot the dependence of the minimum temperature on V g 
for a range of quality factors, which shows that final tem- 
peratures may be increased somewhat by increasing the 
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VI. DISCUSSION AND CONCLUSION 



The results obtained above are consistent with heuris- 
tic arguments. The response of cooling to the measure- 
ment strength is as expected: for very weak continuous 
measurements, we do not learn enough about the state of 
the system to cool it effectively, and can in fact heat the 
system due to acting on our poor information. For very 
strong continuous measurements, we gain sensitivity, but 
inject more quantum back-action, and approach a mini- 
mum only asymptotically. The range of improvement is 
limited, however, and beyond a few volts, the benefits 
may not warrant the additional effort. 

Higher drain-source voltages provide a larger signal- 
to-noise ratio, and therefore improve cooling. However, 
since we do not know exactly how our approximations 
will fail as V^ s approaches e/Cs, and we lack a complete 
theory of the SET once more than 2 island states need 
to be taken into account, we have chosen to stay below 
that limit. 

We have made a few additional simplifying assump- 
tions, as a way to indicate a goal, rather than an im- 
mediately achievable experimental realization. First, we 
have assumed a perfectly efficient (and infinite band- 
width) measurement - that is, that no electron passes 
the detector without being detected. While detection ef- 
ficiency is not as much of a problem here as in optical 
experiments, detectors will necessarily be inefficient to 
some extent. Second, we have assumed perfect, noise- 
less feedback. In reality, the actuating gate applying the 
feedback voltage will not provide a perfect noiseless volt- 
age. Also, we have assumed that the actuating gate does 
not affect the SET. This last assumption is realistic, how- 
ever, for two reasons. First, the resonator itself acts as 
a shield between the gate and the SET. Second, since 
the observer knows the voltage on the feedback gate, she 
can subtract that effect off the SET signal, while adding 
some noise. 

As mentioned previously, a quantum mechanical har- 
monic oscillator and a classical one are indistinguishable 
as long as the wave function is Gaussian, which is the 
case in the present analysis. Therefore, although the os- 
cillator is near the quantum mechanical ground state, the 
SET measurement of position will not show any quantum 
behavior. In the face of these limitations, it is a pleas- 
ant result that experimentally obtainable situations to- 
day allow for the feedback cooling of an resonator to the 
point that quantum behavior could become distinguish- 
able from classical behavior with an appropriate mea- 
surement scheme. 
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APPENDIX A: THE SPECTRA OF THE SET 
SHOT NOISE AND BACK-ACTION 



Here we discuss briefly how the expressions for the shot 
noise and back-action of the position measurement via a 
SET are obtained. For more details the reader is referred 
to Zhang and Blencowe 8 , (from which we obtain most of 
the following expressions) and Korotkov 3 . 

The SET consists of a central island, which electrons 
tunnel in and out of via junctions on either side. If one re- 
quires that the spacing between the energy levels of the 
electron states on the island are sufficiently large com- 
pared to the voltage drop across the SET, then only two 
island states will be appreciably populated, these being 
the states in which there are n and n+1 electrons on the 
island, for some n. This is because the transition rates 
which connect these states to the other states are sup- 
pressed. The value of n can be set by biasing the central 
island. In particular, n is determined by the condition 



n < 



V ds /2) <n+l. 



(Al) 



As a result, we can write a master equation for the prob- 
ability density for the occupation of the two states. De- 
noting this density by a = (a(n),a(n + 1)) T , we have 



da 
~dl 



—a(n) b(n + 1) 
a(n) -b(n + l) 



(A2) 



and 



where a(n) is the transition rate from n to n + 1, 
b(n + 1) is the transition rate from n + 1 to n. 

If we denote the tunneling rates into the island across 
the source junction and the drain junction (see Fig.^) as 
a_(n) and a+(n) respectively (the plus and minus sub- 
scripts record whether the tunneling event has a positive 
or negative contribution to the SET current), and out of 
the island as 6+ and respectively, then 



a(n) — a+(n)+a-(n), 
b(n + l) = b + (n + l) + b-(n 



!)• 



It is also useful to define 



/(») 
g(n + l) 



a + (n) 
b+(nA 



- a_(n), 
1) - 6_(n 



(A3) 
(A4) 



(A5) 
(A6) 
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In what follows we will repress the arguments of these 
functions, so that a = a(n), b = b(n+l) etc. The solution 
to the master equation is 



cr(i) 



b b 

a a 



a 
-a 



,-(a+b)t 



*(0) 
(a + b) 



(A7) 



From this it is straightforward to calculate the average 
steady state current flowing through the SET, the noise 
spectra of the current, Si(u>), along with that of an arbi- 
trary function, 4>(ri), of the island electron number S^{iJ), 
and their mutual correlation spectrum, C(uS). The aver- 
age current is 



I=(ag + bf)/(a + b), 
and the spectra are 

2ab [<p(n) - (j){n+ l)] 2 



S<j,{u) 



(a + b) (a + b) 2 + us 2 



2e 2 C 2 



(a + b)C 2 



ab 



(f-g)(a 2 g~b 2 f) 



(a + 6) 2 



(A8) 



(A9) 



(A10) 



(ag + bf) 2 (a-b) 2 + ^{ag-bff 



Aab [ab [{a + b) 2 + lo 2 } + (f - g){a 2 g - b 2 f)] 

. (AU) 

Since the force from the island on the resonator is given 
b}£ 

F = Cs{ fcld Cs) [c{Vds ~ 2Vs) ~ ne]2 ' (A12) 



using Eq. (|A9|) we have 



abe 2 C 2 g (2C ~ C g ) 2 [C(V ds - 2V S ) - e(2n + 



2(a + b)Cld 2 (a + b) 2 + lo 2 

Recall that in deriving these expressions we require 
that the two-level approximation is valid, and this de- 
mands that 



V ds < e/Cs, 
k B T « eV ds . 

The tunneling rates are given by 

(An±ty/(gjgg) 



(A14) 
(A15) 



a±(n) = 



6±(n+ 1) 



1 — exp 



-(An±V ds )/f 



(-An±V ds )/(R j Cx) 



1-exp -(-An±V ds )/T 



where 



An 

Vds 

f 



CgVg 

c£v ds 

2e ' 
C s k B T 



Co Vds 



2e 



1 

2' 



(A16) 



Note that the condition which determines n (Eq. IjAlfl l 
is equivalent to —0.5 < An < 0.5. 

From the expressions for the noise spectra we see that 
both sources of noise are effectively white (independent of 
uj) so long as ui 2 is much less than [a(n)+b(n+l)] 2 . If this 
is the case then the simple quantum theory of continuous 
position measurement presented in the main body of the 
paper provides a good model for the SET measurement. 
Note that the actual back-action noise on the position 
of the resonator is the force noise filtered through the 
resonator spectral function. This is therefore 



S F (u)/m 2 eS 



( W 2 



,2\2 



,2, .,2 



IQ 2 



(A17) 



and has the same form as that predicted using the quan- 
tum mechanical model (Eq. (|11|0 . so long as the force 
noise is white. 

We must therefore evaluate \a(n) + b(n + l)] 2 for the 
range of parameters of interest, and verify that it is much 
larger than uj 2 over the relevant frequency range. First 
we note that the form of the spectral equations is such 
that they are periodic in the gate voltage. That is, the 
values of a(n) and b(n + 1) depend only on An, not on the 
particular value of n in question. As a result we merely 
need evaluate \a(n) + b(n + l)] 2 for a single value of n, 
and check all values of An between -0.5 and 0.5. 

Substituting in realistic parameter values (those that 
we use in our examples in the body of paper) in Eas. ljA16|) 
and J!]])), we find that, regardless of the value of An, 



[a(n) + b{n + 1)] > 2 x 10 1 



(A18) 



for the range of initial temperatures that we consider, 
and this is much greater than the range of u> relevant for 
the dynamics of the resonator, as required. Thus, we can 
drop u from the expressions for the spectra, Eas. (|A10p . 
(|AT3)l and fKH) . and use these to determine the param- 
eters fc, (3 and k in the model of the quantum position 
measurement. 
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